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Abstract 

The ID flow of a continuous beam of Bose-Einstein condensed atoms in the presence of an obstacle 
is studied as a function of the beam velocity and of the type of perturbing potential (representing 
the interaction of the obstacle with the atoms of the beam). We identify the relevant regimes: 
stationary/time-dependent and superfluid/dissipative; the absence of drag is used as a criterion 
for superfluidity. There exists a critical velocity below which the flow is superfluid. For attractive 
obstacles, we show that this critical velocity can reach the value predicted by Landau's approach. 
For penetrable obstacles, it is shown that superfluidity is recovered at large beam velocity. Finally, 
enormous differences in drag occur when switching from repulsive to attractive potential. 
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1 Introduction 



The rapid progresses in the technology of guiding cold atoms (using hollow optical fibers [jy], mag- 
netic guides H and microchips |J) opens up the prospect of similar studies of guided Bose-Einstein 
condensed beams, i.e., of guided continuous atom lasers. Indeed, important progress in this direc- 
tion is presently being made, which uses the techniques developed for guiding cold atoms: wave 
guides have been designed for a Bose-Einstein condensate using a blue-detuned hollow laser beam 
[BJ|; Bose condensation has been obtained over a microchip |5|, |j and a Bose-Einstein wave-packet 
has been propagated in a microfabricated magnetic wave-guide Also, a continuous beam of 
cold atoms has been loaded into a magnetic guide, as a first step in order to perform evaporative 
cooling and condensation in the guide |§. 



In the present work we address the question of superfluidity of a continuous (guided) atom laser, 
namely: what are the conditions for the flow to be dissipationless ? A criterion for superfluidity has 
been proposed long ago by Landau M which states that dissipation doesn't occur if the velocity 
of the flow is lower than the critical value v cr a = min {E(q)/q}, where E(q) is the energy of an 
excitation with momentum q. 

Many experiments have been done in liquid helium II trying to test Landau's idea, and indeed 
one finds a critical velocity, but, in many instances it is much lower than Landau's expectation. 
As Feynman first suggested flQl , this is linked to vortex formation, i.e., to non- linear perturbation 
of the fluid (and not to elementary excitations as implied by Landau The experiments done 

at M.I.T. for Bose-Einstein condensates confirm this view: in these systems there is also a critical 
velocity [12, 13], it is lower that Landau's expectation, and it is also linked to vortex formation 



In the following we are concerned with one-dimensional flows, which are relevant to atom laser 
physics. We devise an alternative phrasing of Landau's argument based on perturbation theory and 
identify its limit of validity. In the generic, non-perturbative case, we confirm that, for obstacles 
represented by repulsive potentials, dissipation begins at a velocity lower than the one expected 
on the basis of Landau's argument and corresponds to emission of solitons (which are the one- 
dimensional analogues of vortices). However, we show that Landau's critical velocity is always 
reached for attractive potentials. Furthermore, following a previous study []i"5| , above Landau's 
critical velocity we identify a new (numerically stable) regime, stationary and dissipative. In this 
regime, the drag exerted on an obstacle can be computed with little numerical effort. Moreover, in 
this regime, we show that at large velocity the drag exerted on a penetrable obstacle goes to zero, 
i.e., superfluidity is recovered. 

The paper is organized as follows : In Sec. 2 we set up the theoretical framework and notations. 
The natural criterion for the breakdown of superfluidity is the absence of drag, and in this section 
we show precisely how the drag can be computed. In Sec. 3 we determine the different types of flow 
and the corresponding drag for an obstacle represented by an external potential (a weak potential 
in Sec. 3.1, a delta-peak in Sec. 3.2, a square well in Sec. 3.3 and a Gaussian potential in Sec. 3.4). 
Finally we discuss our results in Sec. 4, where we emphasize the important differences between 
attractive and repulsive potentials in the non-linear regime. 
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2 A criterion for superfluidity 



We work in a quasi one-dimensional regime, or more precisely we use an adiabatic approximation 
where the condensate wave-function \P(r, t) can be cast in the form 16, 15 j 



#(r,i) = ^(x,t)(f>(f ± ;n) , (1) 

where ip(x,t) describes the motion along the axis of the laser, (j) is the equilibrium wave function 
(normalized to unity) in the transverse (f±) direction; it depends parametrically on the longitudinal 
density n(x,t) = f d 2 r±\^\ 2 = \ij:(x,t)\ 2 . The beam is confined in the transverse direction by a 
trapping potential V±(f±). The adiabatic approximation assumes that the transverse scale of 
variation of the profile is much smaller than the longitudinal one. The transverse degrees of freedom 
are not completely frozen, but adapt to the smooth longitudinal dynamics: this is the essence of 
the parametric dependence of 4> on n(x,t). This represents a significant improvement to what is 
generally defined as quasi one-dimensional approach and results in a non-typical nonlinearity of the 
ID reduction of the Gross-Pitaevskii equation below (Eqs. (§J3|)). 

In the regime in which Eq. (|l]) holds 17], the equation governing the time evolution of i[)(x,t) 
reads |l5| (we set units such that h = m = 1): 

-\%z1> + { W + eKM)]} = iM. (2) 

In (g), Vu (x) represents the effect of the obstacle. We restrict ourselves to the case of a localized 
perturbation with lim a; _ i .± 00 V(x) = 0. Such an obstacle can be realized by crossing the trajectory 
of the atom laser with a detuned optical laser beam whose waist is large compared with the 
perpendicular extension of the condensed beam. An other possibility is to bend the trajectory of 
the guided atom laser; this results in an attractive effective potential proportional to the square of 
the curvature (see fll5fl ). 

e(n) is a non-linear term describing the mean-field interaction inside the beam, and the way it is 
affected by the transverse confinement. For a transverse confining harmonic potential of pulsation 
uj± one has (see [15[| ) : 

e(n) = 2u± na sc in the low density regime na sc < 1 , . . 

e(ri) = 2u>j_ ^Jna sc in the high density regime na sc 3> 1 , 

where a sc denotes the s-wave scattering length of the two-body inter-atomic potential (we consider 
only the case a sc > 0, i.e. a repulsive effective interaction). In the following we use a formalism 
allowing to treat both the high and the low density regime, since both will be of interest in future 
guided atom laser experiments. 

We want to characterize the superfluidity of the flow past the obstacle. To this end we compute 
the drag exerted by the atom laser on the obstacle: a finite drag implies dissipation, whereas 
Fd = corresponds to a superfluid flow. Fa is defined as 

r +0 ° dVJx) 

F d (t) = J dxn(x,t)—^. (4) 
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This definition is quite natural: the force exerted on the obstacle is the mean value of the 
operator dV\\(x)/dx over the condensate wave- function. It is rigorously justified by the analysis 
below in term of stress tensor (Eq. @). 

For analytical determination of the drag, we use the following procedure: the ID version of the 
stress tensor of the fluid is 1£] 
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T(x,t) = — Im(^j* dtip) + -\d x ^\ 2 — e[n(x,t)] — V\\(x)n(x,t) where e{n) = / e(p) dp , 

(5) 

and its impulsion density is J(x, t) = Im (^*d x i^) (in our units, it is also the current density). By the 
conservation equation dtJ+d x T+n d x V\\ = 0, the total impulsion of the beam P(t) = dx J(x, t) 
is related to by: 

^ =T(-oo,t)-T(+oo,t)-F d (t) . (6) 

The physical content of (|6|) is clear: dP/dt equals the total force exerted over the beam. One 
part of this force (— Fd) is due to the potential, the other one is the stress on the boundaries of 
the beam (at left and right infinity). Hence Eq. (^) confirms the heuristic guess (||); besides, it 
allows to determine the drag in a simple fashion in the stationary regime where T and P are time 
independent: 

Fd = T(— oo) — T(+oo) in the stationary regime . (7) 

Hence, in the following, we devote particular attention to stationary solutions of @. They are 
of the form ip(x,t) = exp{— ipt} A(x) exp{iS(x)}, with A and S real functions. The density is 
n(x) = A 2 (x), the velocity v(x) = dS/dx and the current J(x) = n(x)v{x) is a constant that we 
note Joe,. From Eq. (||), the amplitude A(x) obeys a Schrodinger-like equation 

' ^ A + [vAx) + e[n(x)] + -4^-1 A(x) = pA(x) . (8) 



2dx 2 { 11 v 1 1 K n 2n 2 (x) 



As discussed in [15], the radiation condition requires that solutions of (Jsj) have no wake far 
down-stream: long range perturbations of the beam only occur up-stream. Hence the boundary 
conditions have to be imposed down-stream: because of nonlinearity one cannot disentangle an 
incident and a reflected part in the perturbed up-stream flow. In all the following we consider a 
beam going in the negative x direction, with down-stream boundary conditions n(x — > — oo) = rioo 
and v(x — ► — oo) = — t>oo (with Voo > 0). Then, the chemical potential is fi = + e(noo). In the 
following we refer to as the beam velocity and to Cqo = (n^ de/dn\ n ) 1//2 as the sound velocity 
(the proper denomination should be "sound velocity evaluated at constant density n^). We also 
express the lengths in units of the relaxation length £ = [2 e(n oc )]~ 1 / 2 . 

For stationary flows, the stress tensor (||) reads 

1 f dA\ 2 J 2 

T(x) = - — +W[n(x)]-V\\(x)n(x) where Win) = -e(n) + an + — . (9) 

2 \ dx J ii 2n 

In regions where the spatial variations of Vii(x) are negligible, T{x) is a constant (as easily seen 
from (||)). 
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3 Determination of the drag 



3.1 Perturbative solution 

Let us first evaluate the drag experienced by the obstacle when the effects of the potential V\\ (x) 
on the flow can be treated perturbatively. In this case, by adiabatically branching the potential, 
one can always find a stationary solution with A(x) = y/n^o + 5A(x) having the correct boundary 
condition (i.e., verifying the radiation condition). Introducing the notation k = 2 [v^ — c 2 ^! 1 / 2 one 
finds (see ||) : 



5A(x) 



f+™ V\\(y) exp{-n\x - y\} dy when Voo < , 

(10) 

2 ^f°° V\\ (y) sm{n(x - y)} dy when > . 



The asymptotic evaluation of (|Io|) when x — > ± oo allows to compute the drag through (0) . From 



(0) and (10), T(— oo) and T(+oo) differ only in the supersonic case where 



T(-oo) = W{ noo ) and T(x -> +oo) = T(-oo) + - ( d6A ^ \ + t f_ S A 2 (x), (11) 



2 V dx ) 2 
with (always in the super-sonic regime) 

6A(x) — » Im (e <Ka: V[,(k)) , (12) 

and = c?xexp(— iAex)Vji(x) is the Fourier transform of Vfi(x). One thus obtains 

F d = when < , F rf = -2noo|Vj|(K)| 2 when Uoo > . (13) 

The gross behavior characterized by (|l^) is general: at low velocity the flow is superfluid, 
whereas at high velocity dissipation occurs. This corresponds to Landau's criterion which deter- 
mines a critical velocity below which the flow is dissipationless: v cr u = mm {E(q)/q}, where E(q) 
is the energy of an excitation with momentum q. For our system E(q) is given by the Bogoliubov 
dispersion relation E(q) = q (c 2 ^ +g 2 /4) 1 / 2 (see, e.g., |l9|]) and the Landau critical velocity is then 
Vcrit = Coo- Hence the present perturbative approach is identical to Landau's criterion since both 
give the same value of velocity for the onset of dissipation and have the same physical content: 
excitation of small, non-localized, perturbations is allowed only above v cr a. 



However Landau's criterion, as well as the reasoning leading to Eq. (13) are, by essence, pertur- 
bative. We show below that, as discussed in the introduction, non-linear effects alter these simple 
perturbative views. 



3.2 delta peak potential 

A first hint of the failure of the perturbative approach can be obtained if studying the effect on the 
flow of a delta-potential Vji(x) = X5(x). In that case a stationary solution is obtained by matching 
two free propagation modes of the laser (i.e., solutions of (||) in the absence of potential) at x = 
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with the condition A'(0 + ) — A'(0~) = 2A^4(0). The upstream and downstream stress tensors are 
constant, they are related by T up = T^own + XA(0)[A' (0~) + A'(0 + )]. Besides, at velocity Voo > Cqo, 
the radiation condition imposes that the down-stream (x < 0) density is constant: this gives 
n(0) = rioo and A'(0~) = 0. Hence, for beam velocity Voq larger than Cqo, in the stationary regime 
(g) yields: 

F d = -2 noo X 2 . (14) 

The same formula would have been obtained by using the perturbative result (|l3l). Hence it may 
seem that the perturbative approach is valid for all range of delta potentials and velocities. However, 
one would expect perturbation to fail for strong potentials and near Voo = Coo , since it is meaningful 
only if |(L4.(x)| <C \/n^, i.e., if |A| <C k = 2 — c^J 1 / 2 (see Eq. (|io|)). Indeed, whereas a naive 
perturbative approach (leading to (|l0|)) predicts that a stationary solution exists for all values of 
Voo, an exact solution in presence of a delta-potential can be found only in some precise range of 
velocity and potentials. This has been studied in Ref. and is illustrated on Figure |l[ Roughly 
speaking, stationary regimes only exist for low and high values of v,^ (in the shaded zone of Fig. |l|). 
In the low velocity case, the stationary profile is a trough (or a bump, depending on the sign of 
the potential) localized on the perturbation ]2Q| , such flows are dissipationless. At high velocity 
the stationary profile has a (non-linear) wake extending to infinity in the upstream direction and 
this corresponds to dissipation (by Eq. (0)). In between, the flow is time-dependent and there is 
no reason why, in this case, the drag should be given by (|i~4]). 




-1 -0.5 0.5 1 



Figure 1: The shaded zone of the plane (voo, X) is the domain of existence of stationary solutions occurring 
for a potential V\\(x) — XS(x). The axis are labelled in dimensionless units. The insets represent density 
prohles n(x)/n oc , typical for the different flows (the condensed beam is incident from the right). Each inset 
is located at values of and A typical for the How it displays. The left (right) lower one is a super fluid 
Bow across an attractive (repulsive) potential. The upper one is a dissipative Row. 

We note here an important dissymmetry between repulsive (A > 0) and attractive (A < 0) 
potentials: the domain of superfluid flow extends up to = Cqo for A < 0, i.e., the critical velocity 
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has the same value as the one predicted by Landau's approach. On the other hand, for A > the 
critical velocity is potential-dependent and decreases significantly. This feature is very general (it 
occurs for all the other potentials we have studied), we will comment on it more thoroughly in 
Sec. 3.4 and in the Conclusion. 



3.3 Repulsive square well 



An other case with analytical stationary solutions is the repulsive square well: Vm (x) is zero, except 
for < x < a where it takes the constant and positive value Vq. In this case (0) yields: 



F d (t) = V n(0,t)-n(a,t) 



(15) 



In the stationary subsonic regime, one can show that n(a) = n(0) (see [|15|]) and the drag is zero : 

For > Cqo instead, ( |i~5| ) yields a finite drag. In this case, 
rioo and n(a) in (|l5[ ) can be computed by quadrature, as a 



this characterizes a superfluid flow 
when the flow is stationary, n(0) - 
solution of 



V2 



a — Lq nint 



dA 



[To-W^ + Von] 1 / 2 ' 



(16) 



where nint denotes the nearest integer, Lq is the period of density oscillations inside the well (which 
is again expressible as an integral) and To = W^oo) — VqHoo is the constant value assumed by 
T(x) inside the well. Of course, (li~6| ) is only valid for a hypersonic stationary solution. In the 
regime where no stationary solution exists, the drag is time-dependent and should be computed 
numerically (as done in Section |3.4| below). 



Figure || displays the evolution of the drag as a function of the beam velocity Voo. The drag has 
been computed only in the stationary regimes. In the subsonic stationary regime it is exactly zero 
pC| ]. In the super-sonic stationary regime the exact result (jl5|) (solid line) has been computed in 
two manners: through the numerical solution of (|j~6|), and also by numerical integration of (||). Both 
methods agree within 4 digits. The exact expression is compared in Figure || with the perturbative 
result (|l3| ) (where the Fourier transform of the potential gives here | Vy (k)\ 2 = (2 Vq/k) 2 sin 2 (K(j/2)). 
For large velocities, the perturbative approach becomes more and more accurate. This is expected 
from (10) since, when kg S> 1 (i.e., at large velocities) the perturbative solution is accurate if 
V < k 2 . 

One can remark in Figure [2] that, even in the supersonic regime, the drag happens to be exactly 
zero for some specific values of Voo [21|. As illustrated in the inset, this occurs when the period of 
the density oscillations inside the well is an exact divisor of its width a. In this case, the density 
is unperturbed outside the well. The modification of density inside the well is minute in the case 
represented in the Figure (4%). However, for stronger potentials it can be quite substantial: for 
a = 2 £ and Vq £ 2 = 5, the maximum density inside the well reaches 2 while the density remains 
unperturbed (i.e., equal to n^) outside the well. 
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Figure 2: Drag exerted by a high density beam on a repulsive square well (Vq £ 2 = 0.2. a — 2£), as a 
function of the beam velocity (F^ and i>oo are expressed in dimensionless units). The dashed line is the 
perturbative result (|L3|). The solid line is the exact drag It is evaluated here only in the stationary 

regime. The inset represents the beam density in the supersonic regime, at a value of velocity (indicated by 
the arrow) where the drag is exactly zero. 



3.4 Gaussian potential 

The generality of the above deductions, based on the study of model potentials (a delta peak and 
a square well) can be tested numerically on more realistic potentials. We now consider the case 
V[|(aj) = V exp{-x 2 /a 2 } (with V > 0). 

As in the previous cases, stationary solutions exist only if the beam velocity Voo is not too close 
to the sound velocity Cqq. In the subsonic stationary regime the density is perturbed only in vicinity 
of the potential and the flow is superfluid (by Eq. (^j). In the supersonic stationary regime, the 
density oscillations extend upstream to infinity, and this corresponds to dissipation. At velocities 
where a stationary regime is possible, we have determined the drag indifferently using (|j) or (^) 
(after having solved ® numerically), whereas in the non-stationary case we used (Q) after having 
solved the time dependent equation (g). The results are presented in Fig. |3|. 

The behavior of the drag in the stationary regime confirms what is expected from the results of 
the previous sections. In particular, the critical velocity for the onset of dissipation is lower than 
Cqq. The reason for this is that, in the region of the repulsive obstacle, the density decreases and 
conservation of flux requires that the local fluid velocity increases. As a result, it may happen that 
Landau's criterion is verified (although < Cqo) because the local fluid velocity reaches the sound 
velocity. For slowly varying potentials (i.e., in the regime a 3> £) this argument was put on a firm 
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Figure 3: Drag exerted by a low density beam on a Gaussian potential (Vo £ 2 = 0.2, a — 0.5 £), as a function 
of the beam velocity. The conventions are the same as in Fig. [|. The dashed line is the perturbative result 
(13). The solid line is the exact drag evaluated in the stationary regime. The circles correspond to the 
drag evaluated in the time-dependent regime. The error bars correspond to the extremal values of the time 
dependent Fd(t) (see the text). 



mathematical basis by Hakim [22]. 

From this, one can infer that for attractive potentials instead, the critical velocity should be 
equal to Cqo since for such potentials the density increases in the region of the obstacle and the 
local velocity decreases accordingly |23[| . This has already been observed in the case of a delta- 
peak potential (see Fig. [l]) and also for an attractive square well (see Ref. [15]). We have performed 
numerical checks showing that the same occurs for an attractive Gaussian potential. We will return 
in the final section to the difference between attractive and repulsive potentials. 

Let us now come to the discussion of the time dependent data (the circles of Fig. |3|). They are 
drawn with error bars: this does not correspond to a numerical uncertainty, but reflects the fact 
that the drag depends on time in the non-stationary regime. The "error bars" correspond to the 
extremal values of the time dependent function Fd(t) (see Fig. ||). The initial condition ip(x,t = 0) 
was taken as the stationary superfluid solution in the presence of the well for the sub-sonic case 
v™ lt = 0.5 Cqo to which a Galilean boost exp{i (v 1 ™* — Voo) x} was applied at t = in order to 
reach the desired value of velocity. Accordingly, in the computation, the drag F^(t) starts from 
and after a set-up time reaches a regime where it oscillates around a mean value (see Fig. ||) 
which is represented by the circles in Fig. ||. The oscillations of Fd(t) around its mean value are 
of interest because they reflect the cause of time-dependence of the flow: the numerics indicate 
that whereas the up-stream flow reaches a quasi-stationary pattern, the down-stream density is 
perturbed by solitons, periodically emitted from the obstacle, that propagate in the same direction 
as the flow with a smaller velocity (such a behavior has already been observed by Hakim |22|). As 
illustrated in Figure ||, the drag decreases when a soliton has just been emitted (this was already 
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noted in two-dimensional flows by Winiecki et al. [24|, where vortex pairs are emitted from a 
moving obstacle). 

In the super-sonic regime, it is also important to discuss the discrepancy between the results for 
the drag computed for stationary flows and for time dependent ones (solid line and circles, respec- 
tively, in Fig. ^). This discrepancy by no means implies that the stationary profile is unstable or 
that the asymptotic time-dependent flow is not stationary. On the contrary, numerics indicate that 
time-dependent flows reach an asymptotic stationary regime for velocities at which such a regime 
exists. Moreover the asymptotic density profile is of the expected type (flat down-stream and oscil- 
lating up-stream). The point is that our specific initial condition tp(x,0) does not asymptotically 
lead (when t — * oo) exactly to down-stream density and velocity which have the same value rioo and 
Woo as the initial flow. This is illustrated in Figure ||. On the Figure one sees that the asymptotic 
down-stream density is about 6 % lower than and a simple numerical check shows as well that 
the asymptotic velocity differs from Voo (is is roughly 3 % higher) . This artifact becomes less and 
less important for increasing velocities: since the perturbative approach is more and more accurate, 
it is clear that down-stream modifications of the solution become minor and that the asymptotic 
solution is the expected one. As a result, when increases, the black circles get closer to the solid 
line in Figure [3|. 

4 Discussion and conclusion 

The above presented results illustrate the well known fact that non-linear effects alter the simple 
perturbative (Landau's) approach for the determination of the critical velocity at which dissipation 
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-100 -50 

Figure 5: Density profile in the case Voo — 1.6 Coo, in the 'presence of a Gaussian potential (Vq(, 2 = 0.2, 
a = 0.5 £_). The plot represents the profile (solid line) after evolution of the initial n(x, t = 0) (dashed line) 
during a time t with v^t/^ = 203.6. The density trough initially located at x = has moved to the left 
and is now at x ~ — 100£. It will asymptotically move to left infinity. It is followed by a density slightly 
depressed compared to the initial value . This depressed density is the asymptotic down-stream density. 

occurs (we saw, however, that the perturbative approach could reach a regime of accuracy at large 
velocity). For repulsive potentials the critical beam velocity is smaller than the velocity of sound Cqo 
(which is here the critical velocity from Landau's criterion) because, in the region of the potential, 
the local fluid velocity can reach values higher than the sound velocity. The onset of dissipation 
corresponds to non-stationary flows with a wake asymptotically extending upstream to infinity, and 
down-stream periodic emission of solitons. In fact, another way of explaining the lower stability of 
the dissipationless flow in the presence of a repulsive potential is to remark that, for sub-sonic flow 
over such an obstacle, the density decreases in the region of the potential, allowing easier nucleation 
of solitons. 

On the other hand, we have shown that, for attractive potentials, stationary dissipationless 
solutions exist up to = c^: Bose-Einstein condensates appear to be excellent supports for 
reaching Landau critical velocity, more appropriate than superfluid helium because, in atomic 
vapors, it is simpler to construct obstacles described by purely attractive potentials. 

We also showed, using analytical an numerical examples, that stationary dissipative profiles 
exist in hyper-sonic beams, provided the beam velocity is large enough. From the numerical 
study, these solutions seem stable, and moreover time-dependent flows tend asymptotically to 
such solutions (when they exist). It is interesting to note that dissipation is drastically reduced 
at very high velocity: i.e., superfluidity is recovered. Such an effect should also exist in higher 
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dimensions for penetrable potentials. It can be understood perturbatively: at high velocity, when 
the perturbative approach becomes valid, the relevant wave- vector (denoted k in Sec. 3.1) is large 
and Bogoliubov dispersion relation becomes exactly quadratic. Hence, in this regime one has a 
matter-wave described by the linear Schrodinger equation. The drag in this regime can be shown 
to be proportional to the reflection coefficient which, as well known, decreases at high energy (in 
any dimension). 

Note that, for the sake of clarity, we have always illustrated our conclusions using rather weak 
perturbing potentials. The reasons for this are twofold: (1) for stronger potentials the super-sonic 
stationary regime that we wanted to emphasize is rejected to higher velocity; (2) the numerical 
effort necessary in the study of the non-stationary regime is decreased for weak potentials since the 
domain of time-dependence is reduced. 

Although for stronger potentials the qualitative results remain the same, very interesting new 
quantitative phenomena occur. In particular, enormous differences in drag can occur when switch- 
ing from repulsive to attractive potentials (provided the potential is strong enough). This is illus- 
trated in Fig. |6|. 




1 2 3 



non-statiinaryVegime (V ^ = 1.0 ) 




v Ic 



Figure 6: Drag exerted by a low density beam on a Gaussian potential (Vo £ 2 = ±1.0, cr/£ = 1.0), as a 
function of the beam velocity. The conventions are the same as in Fig. ^. The main figure displays the drag 



for a repulsive potential (solid line) together with the perturbative result (13) (dashed line) which is not 
affected by the sign of the potential. The inset displays a blow-up of the main figure allowing to see the 
(very small) drag for an attractive potential. 

One sees in the Figure that the domain of existence of stationary solutions is markedly different 
between attractive and repulsive potentials and that the values of the drag differ by more than two 
orders of magnitude ! The physical explanation for this phenomenon is subtle. At high velocity, 
contrarily to intuition, the density decreases (increases) in an attractive (repulsive) potential |25f |. 
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Hence, an attractive potential creates a density trough which, being the super-sonic analog of a 
gray soliton, does not create large perturbations in its wake. This results in a small drag. This 
feature is clearly missed by the perturbative approach which gives a drag insensitive to the sign of 
the potential (see Eq. (fb|)). 

From the result of Fig. |6|it would be very interesting to redo with an attractive perturbing poten- 



tial the experiments which have been done at M.I.T. with a repulsive potential (T^, 13]. Although 
the present computations are valid in a quasi one-dimensional regime (whereas the experiments 
were done in truly three dimensional systems) the present results leaves no doubt that the critical 
velocity for the onset of dissipation should increase and that the energy transfer rate from the 
obstacle to the fluid (i.e., F^Voq) should drastically decrease in case of an attractive obstacle [26]. 



Finally, we note that the present discussion sheds some light on the theory of wave resistance. 
The wave resistance is the part of the drag experienced by a body moving in a medium which is 
caused by excitation of waves in this medium (typically surface waves in the case of boats). In 
a superfluid a T = (i.e., in the present work), this is the only source of drag (if one broadens 
its definition in order to include non-linear effects such as vortex or soliton formation). Recent 
experiments of moving spheres in silicone oil [27] have shown that the wave resistance at Landau 
threshold has a smooth behavior as a function of the velocity, contrarily to the expectation of 
perturbation theory | |28| |. The same behavior was found here : non-linear effects smoothen out 
the unphysical step in drag predicted by perturbation theory. This smooth behavior was already 
observed in the experiments done at M.I.T. (l^, [l^] and in numerical simulations by Frisch et al. 
]p9| ] and Winiecki et al. [24] when moving an impenetrable sphere in a superfluid. Hence, Bose 
condensed systems offer an interesting testing ground for ubiquitous non-linear hydro dynamical 
effects, in a particularly simple theoretical framework (the Gross-Pitaevskii equation). 
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